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Abstract 

We use Quantum Monte Carlo methods to determine T = Green functions, 
G(r,u), on lattices up to 16 x 16 for the 2D Hubbard model at U/t = 4. 
For chemical potentials, jx, within the Hubbard gap, < fj> c , and at long 
distances, f, G(f,uj = //) ~ e~l*1/fr with critical behavior: £j ~ |/i — ^ c \~ u , 
v = 0.26±0.05. This result stands in agreement with the assumption of hyper- 
scaling with correlation exponent v = 1/4 and dynamical exponent z = 4. In 
contrast, the generic band insulator as well as the metal-insulator transition 
in the ID Hubbard model are characterized by v = 1/2 and z = 2. 
PACS numbers: 71.27.+a, 71.30.+h, 71.10,+x 
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At zero-temperature a continuous metal-insulator transition driven by a change in chemi- 
cal potential may be characterized by the compressibility, Xc or the Drude weight, D. In the 
Mott insulating phase both D and \ c vanish while they remain finite in the metallic phase 
]T|Q]. In order to describe the metal-insulator transition from the insulator side, we consider 
the zero-temperature Green function G(r,u) 0. At long distances, \r\, and for values of the 
chemical potential, //, within the charge gap, G(r,u = /i) ~ e~^^ 1 . The metal-insulator 
transition may be characterized by the divergence of as the critical chemical potential, /i c , 
is approached from the insulating phase. £j may be interpreted as the localization length 
involved in transferring a particle over a distance r from the electronic system to the heat 
bath lying at energy /i within the charge gap. Under the assumption of hyperscaling, the 
above quantities are expected to satisfy the scaling relations: 

& ~ A~ u , X c ~ A~" {z - d \ D ~ A" {d+z - 2 \ (1) 

where A = — [i c \, d is the dimensionality and v (z) corresponds to the correlation length 
(dynamical) exponent Since the control parameter, A, corresponds to the chemical 
potential, one obtains the additional constraint vz = 1 as well as 5 ~ A"( d+ *) _1 , 5 be- 
ing the doping concentration. In this letter, based on a recently developed numerically 
stable Quantum Monte Carlo (QMC) algorithm to calculate zero-temperature imaginary 
time Green functions [|J, we calculate £j for the two-dimensional repulsive Hubbard model 
at U/t = 4. We obtain: fi c = 0.67 ± 0.015 in units of the hopping matrix element and 
£i ~ |/i — fi c \~ u with v = 0.26 ± 0.05. On the other hand, the compressibility data of Fu- 
rukawa and Imada |J in the metallic phase is consistent with Xc ~ \fJ> ~ ^cl -1 ^ 2 - When 
hyperscaling describes the transition, the above scaling of the compressibility also leads to 
v — 1/4 0. Comparison of those results puts the hyperscaling assumption on a firmer basis. 
The present estimate of v is a more direct determination of the characteristic length scale. 
In contrast, the Mott transition in the one-dimensional Hubbard model satisfies the scaling 
relations (|I|) with v = 1/2 and z = 2 [^]|)|-^]. The generic band insulator in all dimensions 
equally belongs to the universality class v = 1/2 and z = 2 |. The above results point out 
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the anomalous character of the Mott transition in the two-dimensional Hubbard model. 
The Hubbard model we consider reads: 

H-iiN=-t£ cU^ + DX(nj; t - \){n u - \) - H^cU^ (2) 

(i,j),a i i,a 

Here, denotes nearest-neighbors, creates (annihilates) an electron with z- 

component of spin a on site i and n^ a = ct^c^. In this notation half-band filling corresponds 
to fi = 0. We start by considering the non-interacting case, U/t = 0. In Fourier space, the 
single particle energies are given by = — 2t(cos(ka x ) + cos(ka y )), a x , a y being the lattice 
constants. In this letter, the length scale is set by: \a x \ = \a y \ = 1. At zero temperature, an 
insulator-metal transition will occur when fi —>■ fi c = At. For those chemical potentials, the 
zero-temperature Green function J3J at u = fi is given by: 




where N denotes the number of sites of the square lattice and the factor 2 corresponds to the 
summation over the spin degrees of freedom. Numerically, one obtains: G(f, uo = fi) ~ e - '^^ 
with critical behavior: £j ~ |/i — /ici" 1 / 2 Q. One may check that this example of a band 
insulator-metal transition satisfies the above scaling relations (|1|) and is characterized by 
the exponent v = 1/2. At finite values of U/t and half-band filling the antiferromagnetic 
Hartree-Fock approximation equally yields the exponent v = 1/2 for £j. However, this 
approximation does not satisfy the hyperscaling assumption. 

The physical interpretation of is facilitated by considering the single impurity Hamil- 
tonian (Fano- Anderson model) |TU|,|TTJ: 

H = H ^ + ^ b + 7= E (4* + ■ ( 4 ) 

Here, creates an electron in the impurity state at the origin and energy e&. The hybridiza- 
tion between the localized state and the band electrons alters the energy of the impurity 
level to the value: E b = e b + -fi ^£ E . We will assume E b > for all k. When all single 
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particle states of the valence band are filled and the impurity single particle state empty, 
the probability amplitude of transferring a band electron at site r to the impurity state is: 



where the normalization factor is given by a" (-E&) = ( Ej) _ eJ )2 ■ Comparison between 

equations (|]) and (Q) show that the spatial dependence of the two quantities is identical. 
£i may thus be interpreted as the localization length involved in transferring a particle over 
a distance r from the valence band to the heat bath (see equation (|3|)), or to the impurity 
state (see equation (^|)). This definition of £j bears some similarity to that applied in finite 
size scaling studies of Anderson localized states p^2|JT3|| . When the imaginary part of the 
Green function vanishes in the insulating phase, & may be used to study the metal-insulator 
transition. 

To obtain an estimate of the critical exponent v for the Hubbard model, we require an 
accurate determination of the critical chemical potential, \i c as well as of the localization 
length, £/. Both quantities may be obtained from the knowledge of: 

G ' ( '- r) = e(T) — <w e( ~ T) »j — • (6) 

where Cf )(T (r) = e TH c^ a e~ TH . Here |\& ) denotes the ground state of the half-filled (/i = 0) 
Hubbard Hamiltonian (|]). The above quantity may be efficiently calculated with QMC 
methods. Since we are at half-band filling, the sign problem does not occur and we are able 
to consider lattice sizes up to linear dimension L = 16, namely iV = 16 x 16 without any 
serious difficulties. The calculation of imaginary time Green function in the zero-temperature 



auxiliary field QMC algorithm |14|15| was first reported by Deisz et al. ||16|| . However, their 
approach does not incorporate a numerical stabilization scheme and they are thus restricted 
to small values of r (i.e. rt ~ 2.5 ). The authors have developed a numerically stable QMC 
algorithm for the calculation of G a (f,T). The details of the algorithm may be found in 
reference H. All our calculations were performed with periodic boundary conditions. 
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From the knowledge of G a (f, r) on an iV-site lattice, we may obtain an estimate of the 
charge gap. We denote by \^f^) the eigenvector of the Hamiltonian H with eigenvalue E% 
in the iV-particle Hilbert space. With this notation, 

G(f= 0,r) = £G„(f = 0,t)4E l«kJ^ +1 >| 2 exp (-r - <)) . (7) 

cr 

for r > 0. Fig. la, plots G(r = 0, r) for a 16 x 16 lattice at U/t = 4. We may obtain a reliable 
estimate of the charge gap for this lattice size by fitting G(f — 0, r) to the form e _AcT with 
A c = Eq +1 — Eq for large values of r. Fig. lb shows A c as a function of linear lattice size. 
Extrapolation to the thermodynamic limit yields: A c /t = 0.67 ± 0.015. This result stands 
in good agreement with the quoted result of Furukawa and Imada ||: A c /t = 0.58 ± 0.08. 
Since in the notation of equation (g) the Hubbard model satisfies particle-hole symmetry at 
\i — 0, the critical chemical potential is nothing but the charge gap: fi c = A c . 

For values of the chemical potential within the charge gap, < \i c , the zero temperature 
Green function is real and may be obtained through the relation: 

/CO 
drG(f, r)e™. (8) 
-co 

The Green function G(f, r) is computed at half-band filling where the sign problem is not 
present and the statistical uncertainty does not grow exponentially with lattice size. How- 
ever, since we are multiplying the QMC data by the factor e Tfl , the statistical uncertainty 
will grow exponentially with increasing values of r for 77/ > 0. For each lattice size L, we 
have considered the largest distance r = (L/2, L/2). For this distance, G(f— (L/2,L/2),r) 
is plotted in Fig. 2. Due to particle-hole symmetry at /x = 0, G(f = (L/2, L/2), r) = 
—G(f= (L/2, L/2), — r). For the imaginary time integration, (see equation (§)) and values 
of the chemical potential < 0.65t, a cutoff of rt = 10 proved to be sufficient for the 
determination of the Green function (see Fig. 2). G(f= (L/2, L/2), u = fi) as a function of 
lattice size is plotted in Fig. 3 for several values of \i. For lattice sizes ranging from L = 6 to 
L = 16, an exponential decay may be seen within the quoted statistical uncertainty. From 
this data, we obtain an estimate of the localization length £j. With the above determined 
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value of fi c , we plot \fx — /i c | versus (see Fig. 4) on a log-log plot. For all considered chem- 
ical potentials, £j/a < which corresponds to our largest considered distance (L = 16). 
The slopes in Fig. 4 correspond to values of the critical exponent v — 1/2 and v = 1/4. 
The QMC data is consistent with v = 1/4 and seems to rule out the possibility v = 1/2. A 
statistical analysis yields: v = 0.26 ± 0.05. 

The choice r = (L/2,L/2) is very convenient since apart from a sign change between 
L = 4n + 2 and L = 4n lattices the exponential decay of G(f = (L/2, L/2), u = fi) as a 
function of L, is not masked by a non-trivial oscillation. For other choices of f we expect 
to obtain the same results since due to the point group symmetry of the square lattice, all 
/c-points in the Brillouin zone contribute to G(f = (L/2, L/2), r). We could however not 
check this point numerically. 

To confirm the validity of our approach, we consider the Mott transition in the one- 
dimensional Hubbard model at U/t = 4. A similar QMC calculation as described above 
for the two-dimensional case yields a value of the correlation length exponent consistent 
with v — 1/2 (see Fig. 5). We have obtained fi c /t = 0.66 ± 0.015 which is consistent with 



the exact result of Lieb and Wu (TTj: fi c /t = 0.643. Chains of linear length up to L = 24 
were considered. This result stands in agreement with the assumption of hyperscaling with 
exponents v = 1/2 and z = 2 HH-fl. 

In conclusion, we have determined the correlation length exponent from the knowledge 
of £j in the insulating phase of the two-dimensional Hubbard model and obtained: v = 
0.26 ± 0.05. Under the assumption of hyperscaling with exponents v = 1/4 and z — 4, 
this result stands in agreement with compressibility measurements in the metallic phase 
0. We have shown that a similar calculation for the one- dimensional Hubbard model, 
yields results consistent with v — 1/2. Several anomalous aspects of the metal-insulator 
transition are inferred when it is characterized by the universality class, v — 1/4 and z = 4 
0. Based on a single-particle theory, the exponents u — 1/4, z — A are consistent with 
the interpretation of the Mott transition driven by the divergence of the effective mass as 
opposed to the vanishing of the number of charge carriers. This statement is supported by 
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the compressibility in the metallic phase || as well as by the high frequency Hall coefficient 
fT8|j . Another consequence of the exponent v = 1/4, is the behavior of the Drude weight 
in the vicinity of the Mott transition: D ~ 5 2 , 5 being the hole-density. As a by-product, 
we have produced an accurate estimate of the charge gap for the two-dimensional Hubbard 
model at U/t = 4: A c /t = 0.67 ± 0.015. From the technical point of view, we have 
introduced an efficient method to obtain information on the nature of the metal-insulator 
transition by approaching the transition from the insulator side. The most important fact, 
is that for models which show particle-hole symmetry, such as dimerized Hubbard models, 
the here presented method is not plagued by the sign problem and large lattice sizes may 
be considered. 
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Figure captions 

Fig. 1 a) In G(f = 0, r) as a function of r for the half-filled (/i = 0) 2D Hubbard model 
at £//t = 4 on a 16x 16 lattice. The solid line corresponds to a least square fit 
of G(f = 0, r) to the form exp(— A c r) at large values of r. b) A c as a function of 
linear lattice size L. The solid circle at 1/L — corresponds to A c as obtained 
in reference ||. 

Fig. 2 \G(f = (L/2, L/2), r)| as a function of system size and imaginary time r for 
the 2D Hubbard model. 

Fig. 3 ln|G(r = (L/2, L/2), uj = /i)| as a function of distance and chemical potential 
for the 2D Hubbard model. The solid lines correspond to a least square fit of 
\G(r = (L/2, L/2), cu = fj)\ to the form exp(-|r|/6) for L > A. 

Fig. 4 Localization length £j versus \fi — fi c \ for the 2D Hubbard model. The solid 
lines correspond to two values of the correlation length exponent: v — 1/4 and 
v = 1/2. The solid circles correspond to the QMC data. 

Fig. 5 Same as Fig. 4 but for the ID Hubbard model. 
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